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Abstract. The proposed SCORE detector consists 
of a large array of light collecting modules designed 
to sample the Cherenkov light front of extensive air 
showers in order to detect high energy y-rays. A large 
spacing of the detector stations makes it possible 
to cover a huge area with a reasonable effort, thus 
achieving a good sensitivity up to energies of about 
a few 10 PeV. 

In this paper the event reconstruction algorithm 
for SCORE is presented and used to obtain the 
anticipated performance of the detector in terms of 
angular resolution, energy resolution, shower depth 
resolution and y / hadron separation. 

Keywords: shower-front-sampling detector event- 
reconstruction 

I. Introduction 

The SCORE detector is designed to study y-rays 
in the energy regime from a few lOTcV up to some 
10 PcV. The detector consists of a large array of light 
sensitive detectors with a half opening angle of about 
30° and is used to sample the Cherenkov light shower 
front of extensive air showers (EAS). Each detector 
contains a few photomultiplier tubes (PMTs) which are 
equipped with Winston Cones, facing towards the sky. 
With a light sensitive area in the order of 1 ir? in the 
individual detector stations it will be possible to sample 
the shower front up to high distances from the shower 
core, thus making it possible to place the stations about 
100 m apart from each other. The signals from the PMTs 
will be sampled by fast Flash-ADCs, making it possible 
to use not only the photon density but also timing 
features for event reconstruction. 

A comprehensive overview of the physical motivation 
and the detector concept can be found in another con- 
tribution on this conference |8|. 

In this paper a scheme for event reconstruction and re- 
sults of simulations testing this reconstruction algorithm 
are presented. 

II. Simulations 

The simulation chain consists of air shower simu- 
lations in COSIKSA [4J (version 6.75), the detector 
simulation sim_score and a reconstruction framework. 
Showers induced by y-rays and protons in the energy 
range from 10 TeV up to 3 PcV have been simulated, 
and the results are analysed using different detector 
geometries. The two detector alternatives discussed here 
are: 
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Fig. 1 . Photon density distribution with the described fit function for 
an EAS induced by a 1 PeV y-particle. For demonstration purposes 
a 20 m sampHng is assumed. 

1) 100 m spacing and 0.5 entrance window 

2) 200 m spacing and 1.5 entrance window 

In each energy band showers with zenith angles of 0°, 
15° and 20° have been simulated, and all results given 
here are averaged values over the different angles (the 
same number of showers has been simulated for each 
angle). 

The detector simulation includes atmospheric absorp- 
tion of Cherenkov light and the anticipated night sky 
background. The detector modules are modelled with a 
realistic PMT pulse shape function (taken from ISj) and 
a wavelength dependent quantum efficiency (taken from 
||3]). The angular acceptance of the Winston Cones has 
been studied with ray tracing simulations and incorpo- 
rated into the detector simulation. 

All resolutions given in here are la values of the 
difference between MC input and the reconstructed 
quantity. 

III. Measured quantities 

Firstly, the lateral distribution of the photon density 
is recorded. It can be seen that the intensity falls off 
exponentially up to about 120 m from the shower core. 
From 120 m to 400 m the light distribution follows a 
power law, until it changes to an exponential function 
again (see Fig. [TJ: 

{Pi exp((iir) for r < 120m 
Qi^f for 120m<r< 400m (1) 
P2 exp((i2 7') for r > 400 m 

The part up to 120 m can only be used if a spacing 
of 100 m is assumed, otherwise there are not enough 
data points in this range. The data points beyond 400 m 
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Fig. 2. Arrival time distribution with parabolic fit function. Shown 
is the simulation of an EAS of a 2 PeV y-pailicle, recorded by a 
detector array with a spacing of 200 m. 
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Fig. 3. Widtli distribution for the same shower as in Fig. |2] with a 
linear fit from 200 m to 500 m. 



usually have a quite poor signal to noise ratio, making 
the power law part the most important for SCORE. 

If continuity at the connection points (120m and 
400 m) is required in the fits, the number of free param- 
eters can be reduced to four: Q, k, di and d2- Typical 
values for the power law index k are in the range of -2 
to -2.6. 

Secondly, the peak (highest entry in signal histogram) 
and the width of each signal (FWHM) are determined. 
Plotting the peak times versus the distance from the 
shower core one finds that a parabolic function describes 
the data very well (Fig.|2]i. The width of the timing signal 
increases linearly at distances greater than about 200 m, 
while at lower distances the width is roughly constant 
due to the time resolution of the PMTs (Fig. O. 

IV. Shower Core Reconstruction 

The simplest approach for the reconstruction of the 
shower core is to calculate the centre of gravity of 
the measured light intensities (used for example in 
AIROBICC, see |6|). While being robust and fast, this 
method can achieve only mediocre resolutions as it does 
not take into account the non-linearity of the lateral 
photon density function (LDF) shown in Fig. [T] 

A better resolution can be achieved if the known 
LDF is fitted two-dimensionally to the recorded light 
intensities. Using the power law part of the LDF (see 
eq. [TJ, a resolution of better than 5 m can be achieved 
at higher energies. Towards lower energies both methods 
become less accurate when the signal approaches the 
noise level. While at higher energies the spacing of 
the detector stations makes no difference, the smaller 
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Fig. 4. Resolution of shower core reconstruction versus y-ray energy 
for two different detector layouts with 100 m (blue) and 200 m (red) 
spacing (see section and two reconstruction methods: centre of 
gravity method (dashed lines) and LDF fit (solid lines). 
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Fig. 5. Energy resolution versus y-ray energy for two different 
detector layouts with 100 m (blue) and 200 m (red) spacing of stations 
and different sensitive areas per station (see section ITU. 



spacing significantly decreases the threshold for a good 
reconstruction (see Fig. |4]l. 

V. Energy of Primary Particle 

Simulations show that the photon density at any point 
of the Cherenkov shower front is proportional to the 
energy of the primary particle, which makes the energy 
reconstruction a straightforward process. However, at 
small distances from the shower core the photon density 
is also affected by the depth of the shower, while at 
large distances the signal to noise ratio gets worse. 
The best trade-off has been found at a distance of 
220 m, independent of the detector geometry. Taking the 
functional value at 220 m of the fitted LDF (eq. [T]), an 
energy resolution of better than 10% can be achieved 
at high energies. At lower energies, the detector variant 
with 100 m spacing performs significantly better due to 
the four times larger number of data points in the same 
distance interval, providing a better signal to noise ratio 
(see Fig. |5]l. 

VI. Direction of Primary Particle 

It can be seen from Fig. |2] that the peak arrival times 
are distributed according to a parabolic function for a 
vertical shower. For non-vertical showers, this parabola 
is tilted by the angle of the incoming particle. For finding 
the direction of the particle, the superposition of a plane 
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Fig. 6. Angular resolution versus y-ray energy for two different 
detector layouts with 100 m (blue) and 200 m (red) spacing of stations 
and different sensitive areas per station (see section 
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Fig. 7. Shower depth resolution versus y-iay energy for two 
different detector layouts with 100 m (blue) and 200 m (red) spacing 
of stations and different sensitive areas per station (see section llTl. 



and a three-dimensional parabola is fitted to the data 
points: 

fix,y)^E,x + Eyy + A{x^+y^) + c (2) 

Ex and Ey can then be used to calculate zenith and 
azimuth angle. 

Both detector layouts achieve an angular resolution of 
better than 0.1° at higher energies, while the detector 
with the smaller spacing has a better performance at 
lower energies again (Fig. |6l). 

VII. Shower Depth 

Several methods for the calculation of the shower 
depth are discussed in the literature. The TUNKA col- 
laboration has shown that both the steepness of the inner 
part of the LDF and the widths of the timing signals 
(see sec. HIH i can be used for this purpose These 
methods are called intensity method and width method 
from here on. Simulations for SCORE show that also 
the time delay of the peak arrival times, which is fitted 
by a parabola (Fig. |2]i can be used, if an accurate inter- 
station time calibration can be implemented (called tim- 
ing method). Combining the different methods improves 
the depth resolution significantly. 

In order to improve the performance of both the 
width and the timing method, all signals within a certain 
distance interval from the shower core can be summed 
up before determining peak position and signal width. 
As the number of detector stations available in a certain 
distance interval is proportional to the distance to the 
shower core, this procedure makes it possible to sample 
the width and the peak times up to high distances from 
the shower core. The size of the distance intervals is 
a critical parameter of the reconstruction and has to be 
optimized for each detector layout and energy range. 

If the signals are still too noisy for a good determina- 
tion of the width, a lognormal function can be fitted to 
the signals (see [2]). Especially the determination of the 
width of the signal is improved significantly by that. 

Figure [7] shows the depth resolution for the two dif- 
ferent detector layouts. The 100 ni detector variant uses 
all three described methods, while the 200 m detector 
can use only the width and the timing method, as there 



are not enough intensity data points at the inner part of 
the LDF. Both detectors achieve a resolution of better 
than about 30g/cm^ at energies above 300 TeV, but 
the detector with the smaller spacing performs better 
at all energies as it can combine all three methods for 
determining the shower depth. 

VIII. y / HADRON Separation 

The power to distinguish air showers of y-photons and 
hadrons is an important contribution to the sensitivity 
of a y-ray detector. The usual measure to quantify this 
power is the quality factor, which is defined as 

QF = (3) 

where ep is the fraction of particles of the type p which 
survive the applied cuts. Obviously, it is desirable to 
have a high chance to identify a y-photon as signal 
{Ey close to one) and a low chance to falsely identify a 
hadron as signal (ehadron as small as possible), thus a 
high quality factor 

Two approaches for y / hadron separation are inves- 
tigated: 

A. Longitudinal development 

It has been noted in section IVIII that the shower 
depth can be obtained in different ways. While they all 
work well with both y and hadron showers, it seems 
that the timing method and the width method produce 
slightly different results for different particles. This is 
believed to be due to the differences in the longitudinal 
development of showers of photonic or hadronic origin. 
While photonic showers tend to emit most light near 
the shower maximum, the light from hadronic showers 
is emitted from a wider range of altitudes. This makes 
the signals of hadronic showers slightly longer and leads 
to an overestimation of the shower depth for hadronic 
showers in the widths method. 

Comparing the two depth values can therefore be used 
as an indication for the nature of the particle. At higher 
energies (where both methods have good resolutions) 
this can be used for a basic y / hadron separation. Figure 
[8] shows histograms of the ratio of the depth values 
obtained by the two methods for y and hadron showers. 
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Fig. 8. Ratio of the depth value.s obtained by timing method and width 
method for protons (red) and y-rays (blue) with energies between 
500 TeV and 750 TeV. The quality factor for a separation at 1 is 
about 1.2. 



A quality factor of about 1.2 can be reached for energies 
above 500 TeV. 

B. Multivariate Analysis 

In order to find a structure in the data that may im- 
prove the power for y I hadron separation, self-learning 
algorithms can be used. So far, a boosted decision tree 
algorithm with e-boosting (for reference, see e.g. [[TJ) 
has been implemented and first tests look promising. 
Further tests with different input variables, algorithm 
parameters and algorithms are under way. 

IX. Conclusion and Outlook 

The simulations conducted so far have shown that a 
shower front detector like SCORE can achieve angular 
and energy resolutions of the same level as current 
imaging Cherenkov telescopes (better than 0.1° and 
10% respectively). Also, a satisfactory shower depth 
resolution of better than 30g/cm^ has been found. 
Currently, the quality factor for background rejection 
is only at about 1.2, but this has to be considered as 
a first preliminary result. Optimizations of the shower 
depth calculation algorithm are in progress (e.g. better 
fit routines and functions, cuts on bad data points), and it 
is anticipated that a better shower depth resolution will 
have a positive impact on the y / hadron separation. 
Furthermore, other analytic and self-learning algorithms 
for background rejection are being explored currently. 



Of the two detector layouts presented here the one 
with the smaller spacing and the smaller sensitive area 
per detector station has been found to be superior in 
terms of reconstruction power at lower energies. How- 
ever, it uses about 25% more PMTs and four times 
more detector stations. Therefore, a reasonable approach 
might be to combine the two layouts, i.e. use a part of 
the detector section with smaller spacing for detection 
of lower energy particles and employ a high energy 
extension for sensitivity at highest energies. Currently, 
further simulations with different detector layouts and 
spacings are under way. 

The aim of this part of the project is to find the 
optimal detector configuration for the desired application 
and calculate a sensitivity curve based on the simulated 
reconstruction power of this instrument. 
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